TNO  Defensie  en  Veiligheid 


Nederlandse  Orgam  satie  voor 
toegepa  st  -  na  l  uurwet  en  sr  h  a  ppe  hj  k 
o  nde  rz  oek  y  Net  he  r  lands  Organic  atio  n 
for  Applied  Scientific  Research 


TNO-rapport 

TNO-DV2  2005  A221 

De  ontwikkeling  van  een  PBPK  model  voor  VX; 
Stand  van  zaken 
VU 13-8 13  en  207C 


Lange  Kleiweg  137 
Postbus  45 
2280  AA  Rijswijk 

www  tno  nl 

T  015  284  30  00 
F  015  284  39  91 
info-DenV@tno  rcl 


Datum 

februari  2006 

Auleur(s) 

Ing.  1  1.C*  Trap 

K  ii  hr  icering  rapport 

Gngenibriceerd 

Vaslgesield  door 

[)rs.  MAVWalstar.  A po dicker  MPII 

Vustgesteld  d.d. 

1 9  dccember  2005 

I  De/c  rxibncenng  wijzigt  itiel) 

Titel 

Qngembrkeefd 

M  an  age  mentuiUrekse  1 

Ongerubriccerd 

Same  n  vat  ting 

Ongerubriccerd 

Di 

j 

Kappuruefcsi 

Ongcrubriccerd 

Bijtagen 

- 

t 

i'xemplaamummer 

i  • 

1  w 

Oplage 

26 

A  an  tat  pagina's 

21  (inch  bijlage.  escl  RDP&  distnbuiiclijsi) 

Aaitlal  bylagen 

- 

Distribution  Unlimited 


A  lie  rechrtn  voorbehauden  Ntets  utt  dit  rapport  tnag  warden  vermenigvutdigd  en/of  apenbaar  gemaakt  di>or  iniddel 
van  dmk,  fotokopic,  micro  film  of  op  wdke  andere  wijze  dan  ook.  zonder  voorafgaande  schriftclijkc  toesicmming 
van  I  NO 

Indicti  dll  rapport  in  opdrachl  van  he  I  ministehe  van  Defensie  vverd  uilgehraehl.  word*  voor  de  reck  ten  en 
vcrplichtmgen  van  dc  opdraehigever  co  opdraehineinei  verwezen  naar  de  1  Model voorwaafden  voor  Onderzockv  en 
\  Jntwikkdmgsopdraehicif  (MVDT  1 W7)  inssen  dc  minister  van  Defence  en  TNO  indreit  deze  op  de  opdrachl  van 
loepassmg  z»jn  verklaard  dan  we  I  dc  bcireffcndc  ict  zake  i  lessen  partijcn  gtsloten  overeenkomst 


AQ  Foi-n-usfo 


O  2006  FNO 


Managementuittreksel 

TNO  l  Kennis  voor  zaken 

lj^ 

De  ontwikkeling  van  een  PBPK  model  voor 
VX;  Stand  van  zaken 
V013-813  en  207C 


P  rob  lee  ms  tel  ling 

Het  /.enuwgas  VX  is  een  van  de  git'tigste 
middelen  die  er  is.  Samen  met  soman,  een 
ander  zenuwgas,  vormi  bet  een  grote  uitdaging 
om  personen  die  mogelijk  worden  bloot- 
gesield  aan  d£n  van  deze  verbind  ingen  van 
een  effect  ieve  therapie  te  Yoorzien.  Her  is  in 
prwcipe  voor  a  Me  strijdgassen  een  enorme 
uitdaging. 

Om  inzicht  te  krijgen  in  de  gedragingen  van 
VX  in  het  liehaam  zijn  diverse  studies 
uttgevoerd  waarbij  kleine  proefdieren  /ijn 
blootgesteld  aan  dit  middel  (intraveneus  en 
op  de  huid)  waama  de  toxicokinetiek  werd 
bepaald.  De/e  methodiek  leven  goede 


in/ichten  op,  maar  is  hewerkelijk  dtxvr  onder 
andere  het  gehruik  van  pmcfriieren. 
Computersimulaiies  van  dergelijke  studies 
kunnen  inzicht  verse  ha  If  en  zonder  dergelijke 
inspannende  experiment  en.  Vooralsnog  zijn 
dergelijke  studies  nog  steeds  noodzakelijk 
vanwege  de  input  en  valid  at  ie. 

In  hei  kader  van  het  programme  V013, 
‘Passieve  verdedtging  tegen  NBC-wapens’ 
is  opgenomen  het  ontwikkelen  en  opzetten 
van  Fysiologisch  gebaseerde  Pharmaeo- 
kinetische  mode  Hen  (PBPK)  we  ike  de 
toxicokinetiek  van  soman  (opgeleverd)  en 
VX  kunnen  besehrijven.  De  ontwikkeling 
van  dergelijke  mode  Hen  staat  met  op  zich 
maar  kan  good  gebruikt  worden  om  onder 
meer  expositiesoenario's  en  therapeutisehe 
regimes  te  beoordelen. 

Beschrijving  van  de 
werkzaamhcden 

Naast  het  opzetten  van  een  computernxxld 
bestaal  een  groot  gedeelte  van  de 
werkzaamheden  uit  het  verzamelen  van  tie 
noodzakelijke  inputparameters.  Naast  open 
bronnen  a  is  de  lileratuur  en  rapporten  is  een 
ilinke  inspan ning  gelcverd  om  via  (bio)- 
chemiseh  ondcr/oek  diverse  parameters  te 
meten.  Van  vorschillende  weefselhomo- 
genaten  van  een  cavia  zijn  experimenten 
uitgevoerd  om  de  hydmlyiisehe  aetiviteit  te 
kunnen  kwamilkeren  en  de  verdeling  van 
VX  tussen  bloed  en  wee  f sets  te  bepalen.  De 


hindmgssnelheidsconstante  van  VX  aan  de 
(voor  de  mens  en  died  essentiele  enzyme n 
AeetyJCholine  esterase  (AChE), 

B lit yrylCho line  esterase  (BuChE)  en 
Carboxyl  esterase  (Cal;)  zijn  gedeellelijk 
bepaald. 

Resullalen  en  conciliates 

Een  gedeelte  van  de  bemxligde  input  - 
parameters  is  vergaard  of  kan  worden 
gesehat.  De  voorspellingen  van  het  model 
met  betrekking  lot  de  toxicokinetiek  van 
VX  zijn  nog  niei  'klopperuT  met  de 
waarden  vergaard  uit  dierexperimente n .  Dit 
betekem  onder  andere  dal  wellieht  nog 
verde r  gesleuteld  moet  worden  aan  hot 
model  'an  sich'  en  dat  de  pogingen  om  de 
inputparameters  te  verzamelen 
gecontinyeerd  moeten  worden. 

Toepasbaarheid 

Met  een  PBPK-model  kunnen  velerlei 
scenario’s  worden  doorgerekend  zonder 
kostbare  dierexperimente  n  Tevens  vomu 
een  dergelijk  rrnxiel  een  goede  basis  om 
uitspraken  te  doen  over  de  effecten  van  een 
intoxicatie  in  de  mens,  En  om  in  de 
loekomst  therapiesirategieCn  re  modelleren. 
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Samenvatting 


In  vervolg  up  de  onlwikkeling  van  het  PBPK  model  voor  soman  in  de  eavia  en  do  mens 
word!  nu  ecn  model  voor  VX  upgesteld.  Pit  rapport  beschrijli  heknopt  de  werkzaamheden 
en  de  stand  van  zaken  om  lot  een  PBPK  model  voor  VX  te  komcn.  Ecn  PBPK  model 
kan  met  dc  fysiologic  als  basis  (Physiologically  Based}  de  lotgevallcn  (Pharmaco  Kinetics) 
van  een  verbinding  na  blootstelling  voorspellen  in  een  dier  of  uiteindelijk  in  de  mens. 

De  toxicokinetiek  van  hel  zenuwgas  VX  is  wczenlijk  anders  dan  die  van  de  zenuwgassen 
sarin  of  soman,  Vooral  het  persisicnte  karakter  van  VX  bij  vergclijkhare  toxische  dosis 
komt  hierbij  naar  voren.  Bij  die  rex  peri  me  men  is  gehleken  dat  VX  nog  zeer  king  na 
blootstelling  in  hel  lichaam  in  loxicologisch  relevante  concentraties  eireuleeru 
Net  als  bij  hel  soman  model  moet  voor  dit  model  een  Oink  aantal  parameters  worden 
vergaard,  hetzij  via  de  lileratuur  danwel  via  eigen  onderzoek,  Als  basis  is  hel  bestaande 
6-compartimenten  model  voor  soman  gertomen,  aangepasi  en  uitgebreid  met  ecn  extra 
eompartiment,  namelijk  de  huid.  In  dil  model  zijn  tevens  de  readies  met  de  relevante 
enzymen:  ACHE,  BuChEen  CaE  af/onderlijk  beschreven. 

Met  de  introductie  van  een  huid  compartiment  ligt  nu  de  mogclijkheid  open  om  dermalc 
opnames  te  simulercn,  we  Ike  met  name  voor  VX  van  be  king  zijn. 

Gebruikmakend  van  de  bestaande  set  aan  input  para  meters  is  het  met  mogelijk  om  de 
toxicokinetiek  van  VX  na  tntravencuzc  toediening  in  de  (behaarde)  eavia  goed  ic 
voorspellen.  Zodra  het  aantal  lang/a  me  bindingsplaalsen  (CaE)  wordt  verhoogd  is  het 
mogelijk  om  dc  toxicokinetiek  van  VX  in  de  haarloze  eavia  goed  te  simulcrcn. 
Voorlopige  resuhaten  van  pereutane  blootstelling  geven  een  redelijk  Mil'  van  de 
toxicokinetiek  van  VX  in  hel  hfoed  in  de  haarlo/.c  eavia  wcer. 
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Summary 


As  a  follow-up  on  the  development  of  the  FBPK  model  for  soman  in  the  guinea  pig  and 
man  a  new  model  is  being  developed  for  VX*  This  report  describes  in  short  the  work 
and  results  which  has  led  to  this  model.  Based  on  the  Physiology  (Physiologically- 
Based)  of  the  investigated  species  the  date’  or  kinetics  (Pharmaco  Kinetics)  of  a  compound 
can  be  predicted  after  an  intoxication  in  several  organs  or  compartments  by  using  a 
computer  simulation  model 

The  toxicokinetics  of  the  nerve  agent  VX  are  fundamentally  different  from  those  of  the 
nerve  agents  sarin  or  soman.  Especially  the  persistent  character  of  VX  in  blood  at  equitoxie 
doses  is  apparent.  Animal  experiments  have  shown  the  presence  of  lexicologically  relevant 
concentrations  of  VX  in  blood  a  long  time  after  intoxication, 

A  number  of  input  parameters  need  to  be  obtained  cither  through  literature  searches  or 
by  research.  The  6-compartment  model  for  soman  was  adapted  and  expanded  with  an 
additional  compartment,  the  skin.  Also  the  reactions  with  the  enzymes,  AChE,  BuChK 
and  CaE  were  described  individually. 

With  the  implementation  of  a  skin  compartment  it  becomes  possible  to  simulate  dermal 
exposure  scenarios. 

Using  this  model  and  the  available  set  of  inpuL  parameters  it  is  not  yet  possible  to 
simulate  the  intravenous  kinetics  of  VX  in  the  "normal  guinea  pig.  After  adapting  the 
numbers  of  slow  binding  sites  (CaE)  the  toxicokinetics  in  the  hairless  guinea  pig  arc 
simulated  very  well.  Some  preliminary  simulations  of  percutaneuos  exposure  also  show 
a  reasonable  111  of  the  VX  toxicokinetics  in  the  blood  of  the  hairless  guinea  pig. 
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I  Inleiding 


In  vervolg  op  de  onlwikkeling  van  een  PBPK*  model  voor  soman  in  dc  eavia  en  tic 
mens  ||  1 1  Trap  2005]  worth  nu  ecn  poging  ondernomen  urn  een  dcrgelijk  model  op  le 
zctten  voor  VX.  Mel  behulp  van  een  dcrgelijk  model  kunnen  de  lotgevallcn  oftewel  de 
loxieokinetiek  van  een  agens  in  hel  lichaam  worden  voorspeld, 

VX  vorml  door  de  hugere  toxiciteil  en  persistenlie  ten  op/iehte  van  andere  'klassickc' 
zcnuwgassen  |  [ 2 1  Van  der  Schans  2003 1  in  het  lichaam  ecn  I  linke  uildaging  om  le 
komen  tot  ecn  effeclieve  Iherapie  danwel  de  onlwikkeling  van  een  profylaclieum,  Hen 
eomputersimulatiemodel  dal  de  loxieokinetiek  kan  voorspellen  en  inziehl  kan  verschalTen 
in  de  sccundaire  readies  zoals  de  hydrolyse  en  de  readies  met  belangrijke  enzymen  kan 
hiertoeeen  belangrijke  bijdrage  leveren. 

Toi  op  heden  zijn  nog  geen  aanwijzingen  le  vinden  voor  een  werkend  model  dal  de 
loxieokinetiek  van  VX  kan  voorspellen.  Niet  allcen  hel  ontbreken  van  voldoende  validate- 
gegevens  zoals  de  loxieokinetiek  in  diverse  relevanie  species,  zoals  de  eavia  en  de 
marmoset,  maar  ook  van  de  in  pul  parameters  die  nodig  zijn  om  ccn  model  le  voeden, 
zijn  hier  debet  aan.  ReeentcJi  jk  is  wel  door  Sweeney  cn  co-auteurs  1 3]  een  model 
bcschrcven  voorde  loxieokinetiek  van  soman.  Dit  vorml  samen  met  het  reeeni  door 
TNG  ontwikkelde  model  voor  soman  1 1 1  een  aardige  basis  om  op  voorl  le  bon  wen. 

De  loxieokinetiek  van  VX  onderscheidl  zich  van  die  van  de  G-sioffcn  zoals  soman  en 
sarin  [Van  der  Schans,  2003J.  Niei  allcen  is  zoals  eerder  genoemd  VX  nicer  persislcnl 
maar  dc  eliminatie  van  dc  isomeren  blijkl  vooralsnog  weinig  stcrcoscleclief  ic  zijn. 
Onderzoehl  zal  moeien  worden  in  hoeverre  dit  valt  le  verklaren  uit  dc  binding  aan 
AChE,  BnChE  en  CaE  en  dc  hydrolyse  van  VX  in  de  wfeetse1s.  A  Is  onderdeel  van  dil 
project  wordi  een  aantal  van  deze  toxicanl-specitlcke  gegevens  verzameld.  Hel 
bestaande  TNO-model  voor  soman  zal  worden  aangepast,  t  mid  at  in  dal  model 
gesproken  wordi  over  dc  binding  van  allcen  de  toxische  P(-)  isomeren  van  soman  met 
AChE  cn  de  eomhinatie  van  BuChE  en  CaE  als  1  enzym  (  Lumped  E').  Het  model  zal 
derhalve  worden  omgebouwd  om  zowel  de  separate  enzymen  le  kunnen  volgen  als  ook 
een  slereoselectievc  hydrolyse. 

Ecn  belangrijke  nieuwe  loevoeging  in  het  model  is  dc  mogclijkheid  voor  perculane 
toediening  van  VX.  Damp-  cn  vloeislofblootsiellingen  vormcn  nJ.  de  grootstc  risico's 
bij  exposilie  aan  VX.  Tevens  is  uit  de  dierex  peri  men  ten  duidelijk  geworden  dal  relatief 
grote  variaties  in  dc  bloedconcentraties  van  dc  VX-isomercn  in  dc  dicrcn  ondcrling 
opireden  na  perculane  bloolstelling.  Ecn  mogelijke  verklaring  hiervoor  zou  gezoehl 
kunnen  worden  bij  de  fysiu logic  cn  dc  biochemisehe  aetiviteit  van  de  hu id. 
Bovenstaande  onderschrijfl  dit,  daamaast  is  er  de  wens  om  naast  de  intravencuzc 
toediening  cn  dc  respiratoire  bloolstelling  een  apart  huidcompartimem  op  te  nemcn  in 
hel  model  om  nieuwe  deze  relevanie  blootstellingsrouie  te  incorporeren  in  het  PBPK- 
model. 


PHPK  is  du  engelse  afkorting  voor  Thy  srologically- Based  PharmucoKtneiic. 
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2  Methodiek 


2.1  Isolatie  van  de  VX  isomeren 

Voor  het  volgen  en  besludcren  van  de  htochemische  aetiviteit  van  de  aparle  isomeren  is 
op  kleine  sc  haul  een  syniheseprogramma  opgestart. 

Circa  6,4  gram  (+)VX  met  cen  optische  /.uiverheid  van  97,9%  (HPLC  gcgcvens)  en  cen 
chcmlsche  /uiverheid  van  97%  (NMR  gcgcvens)  werd  vcrkregen.  De  andcre  isumeer  (-) 
VX  werd  met  cen  ophrengst  van  circa  8,0  gesynthctisecrd  met  een  optische  /uiverheid 
van  98,4%  (HPLC  gegevens)  en  cen  ehemische  /uiverheid  van  98,5%  (NMR  gegevens). 

2.2  Bepalen  van  de  hydrcdysesnclheid  van  de  isonieren 

De  enzymaiische  hydrolyse  van  VX  in  wcef'selhomogenaat  afkomstig  van  under  andcre 
caviars,  marmoset,  rhesusaap  en  een  stukje  lever  van  de  mens  werd  bepaald  door  de 
conccntratie  VX  te  meten  na  incuhatie  met  het  beireffendc  weefselhomogcnaai.  Oni  le 
onder/oeken  of  deze  hydrolyse  stereoselectief  plaatsvindl  warden  de  homogenaten 
gei'neubeerd  met  de  gcTsoleerde  isomeren  van  VX  en  gevolgd  in  de  lijd. 

Hen  serie  van  circa  15  Eppertdorf  vaaijes  (2  ml)  werd  gevuld  met  met  cen  portie  wee  f sc  I - 
homogenaal  (500  u!  in  25%  in  saline)  en  werd  gemengd  met  (+)-VX  of  (->VX 
(eindconcentratie  10  ug/ml).  De  Eppcndorf  vaatjes  werden  in  een  waterbad  geincubeerd 
bij  57  °C.  Op  bepaalde  tijden  werd  cen  vaalje  uitgenomen,  de  interne  standaard 
(n-propyl-VX,  50  pi  a  lOOug/m!)  toegevoegd,  20  ui  IN  NaOH  toegevoegd,  en  word  het 
VX  geextraheerd  met  tweemaal  I  ml  hexaan.  Vervolgens  werden  de  hexaanlagen  samen- 
gevoegd  en  gcanalysecrd  met  GC-NPD. 

De  gaschromalograaf  (Fisons,  Milaan  Italic)  was  uitgerusl  met  een  NP  detector  en  een 
split-split  less  injector  Tempera  tuur  van  de  injector  en  detector  waren  260  en  300  LJC, 
respecticveiijk.  Tijdcns  de  injeclie  werd  de  splilleiding  voor  60  s  gesloten.  De  stationaire 
I'ase  van  de  kolom  was  cen  CPSiL  8  (I  =  25  m,  i.d.=  Ot25  mm.  Him  =  1,2  urn.  Varian, 
Middelburg).  Het  ovenprogramma  was  als  volgt  ingesteld:  100  <fC  gedurende2  min, 
gevolgd  door  een  stijging  met  25  uC/min  naar  250  °C,  wclke  werd  aangehouden  voor 
5  min,  gevolgd  door  een  stijging  met  25  °C/min  naar  300  "C, 

2.3  Bindmgssnclhcidsconstaiiten  met  AChE,  BuChE  cn  CaE 

Voor  de  bepaling  van  de  bmdingssnelheidsconstanten  (k)  van  dc  VX  isomeren  met 
AChE  en  BuChE  werd  gebruik  gemaakt  van  eavia  AChE  en  humaan  BuChE.  Eerder 
werd  door  Langenberg  en  De  Vette  de  bindingssnelheid  van  VX  met  CaE  vaslgestcld. 
Humaan  BuChE  was  verk regen  van  een  externe  bron  en  eavia  AChE  is  getsoleerd  nil 
cryrh  rc  >ey  ten .  Een  maa  I  geiso lee  rd  we  rd  een  c  n  zy  mat  i  sc  he  k  I  e  u  ri  n  gsreaet  ie  (eo lo ri  met  ri  sc  he 
reactie,  de  z.ogcnaamde  Ellman  reaetie)  uitgevoerd  met  de  VX  isomeren  om  dc  k-waarde 
Le  bepalen. 

Reactievergelijkingen  van  de  enzymkineiiek 


E  +  I  — >  |E  l|  — >  El 


(I) 
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ACME  of  BuChE 

Aeelyllhiocholinejodide  (ASChJ)  — »  Thiocholinc  +  Acclaat 


(2) 


BuChE 

Butyryllhfocholinejodid  ( BuSChJ)  — *  Thiocholine  +  Acclaat  (3) 

Thiocholinc  +  DTNB  — »gdc  klcur  (2-nitro-S  thio-bcnzoaat)  (4) 

Deze  gele  kleur  heel!  een  absorptiemaximum  bij  412  nm  cn  kan  met  de  LJ  V/Vis 
spectrometer  worden  bepaald. 

2.4  BloedAveefsel  verdeJing  van  VX 

Weefselhomogenalen  (25%  w/w)  worden  gemaakt  met  een  veronalhulTer  (niei  nadcr 
beschreven)  vanwegc  de  hydrolytische  albraak  van  VX  door  andere  buffers.  Hot  is 
cssentieel  oni  de  homogenaatoplossing  up  pH  7,40  (±0,05)  le  hrengen  in  verhand  met 
de  pKa  van  circa  9  van  VX.  Gnder  fysiologtschc  omslandigheden  is  VX  gedeeltelijk 
ge laden  helgeen  de  part  it  ie  bepaling  flink  bcmoeilijkt. 

De  bind  ingspl  nut  sen  in  de  wee  f sets  worden  bezel  door  een  grote  overmaat  aan  remmer 
toe  te  voegen,  vanwege  de  hoge  alfinileil  van  soman  voor  de  AChE,  BuChE  cn  CaE 
wordt  deze  gebruikt.  EDTA  wordi  loegevoegd  om  de  enzymatische  hydrolytische 
aetiviieit  te  vcnmmdercn". 

Bovenslaande  acties  worden  uitgevoerd  om  lal!e‘  mctabolc  aeiiviteil  van  het  homogenaat- 
mengsel  plat  te  leggen.  Vervolgens  worden  aan  de/e  homogen aaloplossingcn  een  kleinc 
hoeveetheid  VX  loegevoegd.  VX  zal  zich  gaan  verdclen  tussen  diverse  lagen  tn  dit 
mengsch  na  even wie his inslel ling  wordt  met  hexaan  de  niel  opgenomcn  hoeveelheid  VX 
gecxtralieerd  en  vervolgens  verder  geanalyseerd  met  GC-NPD,  Door  de  resultalen  op 
clkaar  le  delen  en  rekening  houdend  met  het  effect  van  de  buffcroptossing  kan  uiteindelijk 
de  parti  tie  of  verdcling  worden  berekend  [  I J  en  [8|. 


Vooralsnog  blijft  de  vraag  of  soman  en  VX  de/ elide  bindingsplaatsen  bc/eilen. 
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3  (Bio-)  Chemische  input  parameters 


3*1  Algemene  parameters 

Ecn  gcdeclle  van  de  inputparameters  zuals  deze  voor  het  PBPK  model  voor  soman  zijn 
gebruikt  kunnen  onverkort  worden  overgenomen  of  daarvan  worden  afgelcid. 


Tabel  I  Algemene  inputparameters  voo r  een  PBPK  mode!  in  de  cavia  (5(H)  g). 


Compartimenten 

Parameter 

A  rte  heel 

Veneus 

Longen 

Arm- 

Rijk- 

Lever 

Huid 

Bloed 

Bloed 

doorbloedpp 

doorbtoed^ 

Volume  fracti e 

0.02 

0T04 

0.012 

0t53 

0.17 

0t06 

0,17 

Volume  (ml) 

10 

20 

6 

275 

85 

30 

85 

Bloed  How  Jractie  # 

1 

1 

1 

0.19 

0,46 

0,30 

0.05 

Alveolaire  flow  (ml/s) 

2 

AChE  cone.  (nM) 

9.7 

4,5 

1.6 

49 

6f3 

1 

BuChE  (nMf1 

21 

21 

11 

0 

4 

79 

0 

CaE  (nM)131 

440 

440 

170 

5# 

170 

20.000 

5 

De  waarden  voor  volume  en  bleed  fr  aches  hei  ‘poorly  perfused'  of  armdoorbloed  (pp)  eomparlimeni  en  die  van  'richly  perfused’  af 
rijkdoorbloed  compart  intent  Crp)  /ijn  result  an  ten  van  gedefinieerde  compart  i  men  ten. )  3 )  Gegevens  geexirapoleerd  uil  Sweeney  (2(H>5).  |5| 
Data  Langenberg  (2<X)5). "  De  aanlallen  bintiingsplaaisen  of  concent raties  CaE  in  het  armdoorbloed  comport imeni  komen  met  overeen  met 
die  /oals  door  Sweeney  wordt  gepubliceerd,  echter  zij  komen  mecr  overeen  met  het  soman  model  De  grijze  cellen  in  de/e  label  /ijn  sehattingcn 
en  moeten  wcllicht  nog  nader  worden  bepaald 


In  PBPK  model len  wordt  algemcen  gepostuleerd  dal  het  lichaam  voor  76%  uil  ecu 
armidoorbloed  gcdeclle  en  24%  uil  ecn  rijkdoorbloed  gcdeclle  hestaat.  Zodra  er  spccifiekc 
compartimenten  worden  gedefinieerd  wordt  dit  verrekend,  zie  ook  bovenstaande  label 
De  huid  samen  met  het  bloed  en  hei  armdoorbloed  comparti mem  moeten  samen  circa 
76%  vormen. 

3*2  Hydrolyse 

Met  de  losse  isomeren  is  in  de  diverse  weefselhomogenalcn  de  hydrolyse 
snelheidsconslante  bepaatd,  zie  tabel  2. 


Tabel  2  Hydrolyse-snelheidsconstanlen  voor  de  VX  isomeren  bepaald  in  weefselhomogenalcn. 


Compartimenten 

Parameter 

A  rte  heel 

Veneus 

Longen 

Arm- 

Rijk- 

Lever 

Huid 

Bloed 

Bloed 

doorbloed^ 

doorbloed'1* 

ktp  (s  ’) 

4,36*1 04 

4,36*  104 

2,61  *10e 

1,00*1  O'10 

3,29*  1 0  s 

4.65*105 

1,00*10'° 

ktm  (s') 

1,13*1 0  s 

1,13*10'5 

2,04*10  B 

1,00*10 10 

3,33*10  5 

4,06*10  B 

1.00*10“ 

- 

kip  is  de  hydrolyse  snelheidsconslante  voor  (+)  VX  en  kfm  voor  die  van  H  VX.  De  gripe  cellen  in  deze  label  zijn  sc  halt  ingen  en  moeten 
wcllicht  nog  nader  worden  hepaald. 


De  Tnetingen  in  dc  homogenaten  van  de  huid,  rest  van  het  lichaam  en  in  het  vet  laten 
geen  hydrolytische  activiteit  /iert.  Voorzichlighalve  is  de  snelheid  op  ecn  erg  Usge 
waarde  gezel,  dit  belreft  de  spontane  hydrolyse. 
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3.3  Rindiiigssncihetd  van  dc  VX-isomeren  aan  AChE,  BuChE  cn  CaE 

Het  bcpalcn  van  doze  bindingssnclhcidseonstanten  hlcck  moeilijkcr  ic  /ijn  dan  werd 
verwachl.  Om  de  bindingsconstanic  van  de  VX  isomcren  me i  cavia  AChE  te  bcpalcn 
werden  ccrst  uit  eaviahloed  dc  erythroeyten  geYsoleerd  en  opgewcrkt,  waama  vervolgcns 
dc  aetiviteit  werd  bepaald.  Van  ccn  cxlcrne  bron  werd  ccn  kleine  hocveclheid  humaan 
BuChE  verkregen  waarmee  vervolgcns  dc  bmdingsexperimenten  konden  worden 
mtgevoerd.  Humaan  BuChE  heefl  ccn  redeiijkc  stereoselecdv  licit  voor  de  binding  van 
VX,  respecticvclijk  6,44*  10s  M'Vs  1  voor  (+)  VX  en  5,6I*106  M‘Vs'J  voor  {-)  VX-  Er 
was  hclaas  geen  cavia  BuChE  voorhanden  om  tc  controleren  of  deze  siereoselcclmicit 
mogelijk  hetzelfde  zou  zijn,  Dc  schatting  zoals  dcz.c  nu  wordi  gemaaki  is  voor  {+)  VX 
ccn  bindingsconstanic  van  5*I05  M  Vs  1  cn  die  voor(-)  VX  5*10h  M  Vs  V 
Door  Langenberg  cn  De  Velle  is  dc  bindingssnclhcidsconstante  van  VX  met  cavia  CaE 
gemeten  cn  bepaald  op  1 ,5*10^  M’Vs  V 


Tabcl  5 

Bindi  ngssnel  heidsconstanie  ( 

±)  VX  me!  AChE.  BuChE  en  CaE  [M  . 

s'). 

AChE 

BuChE 

CaE 

<+)  vx 

5,39*  104 

5*1 0s 

1.5*103 

(■)  VX 

7,97*1 06 

5*10® 

1.5*103 

Dc  waarden  voor  dc  binding  met  cavia  UuChE  tgrijze  gehied)  /  ijn  sehallingcn  cn  gcbaseml  op  dc 
verhouding  humaan  BuChE  en  ‘Electric  Eel*  AChE, 


3,4  Bloed/weefsel  verdeling  van  VX 

Het  blijft  mocilijk.  om  de  verdeling  van  intact  VX  lusscn  blocd  en  lucht  cn  tussen 
wccfsel(homogenalcn)  en  bloed  tc  bcpalcn.  Niel  alleen  is  VX  rcaclief  (bindingsreacties 
en  hydrolyse)  maar  kan  de  lading  die  het  sukstofatoom  krijgl  in  walcrig  milieu  dc 
‘recovery1  llink  parten  spelcn  waardoorcen  en  ander  mocilijk  interpreteerbaar  wordi. 
Om  toch  het  model  tc  kunnen  voeden  wordi  in  eerste  i  n  slant  ic  gewerkt  mcl  ecu  aantal 
schattingcn  van  deze  parameters.  Momcnteel  worden  de  laatste  wcrkzaamheden  verrichl, 
dc  resultaten  zullen  zosncl  mogelijk  in  het  model  worden  opgcnomen. 
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4  Het  PBPK  model 


4.1  Stmctuur  van  het  model 

Er  wordt  (voorlopig*)  gebruik  gemaakt  van  de  programmeertaal  "Berkeley  Madonna", 
afkomstig  van  Macey  en  Oyster  (University  of  California  at  Berkeley,  CA,  USA), 
waarmee  ook  ervaring  is  opgedaan  voor  het  programmeren  voor  het  soman  model.  In 
figuur  1 ,  is  de  stmctuur  van  het  model  schematisch  weergegeven.  De  huid  is  nu  als  apart 
compart iment  weergegeven,  tevens  zijn  er  nu  drie  verschillende  toedieningroutes 
beschikbaar* 


Hguur  I  Schematische  weergavo  van  hot  PBPK  model  voor  VX, 


Er  is  niet  gekozen  om  de  huid  verder  te  compartimenteren,  zie  ook  figuur  2  voor  een 
representatie  van  een  verdeling  over  de  huidiagen,  Onder  andere  door  het  ontbreken  van 
valide  (input)  gegevens  van  de  huidlagen  is  het  vooralsnog  niet  zinvol  deze  verdeling 
verder  door  te  voeren. 


Mol  hoi  worken  met  Berkeley  Madonna  is  veol  ervaring  opgedaan  do  mogelijkhcden  voor 
uiihreidingen  zijn  beperki,  met  name  op  het  vlak  van  multi-scenario  analyses  en  parameter 
optimal!  satie.  Een  mogelijkhoid  /ou  kunnen  rijn  overstappen  naar  een  ander  pialvotm.  waardoor  hot 
programma  ACSL  eventueel  woer  in  aanmerking  kan  komen, 
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Figuur  2  Sehcmatische  wcergave  van  eon  opname  in  do  huid  on  do  diverse  hmdtagen. 

Er  is  we  I  een  aantal  gegevens  bekend  over  de  huid,  zo  is  volgens  Bookout  [6]  de 
onderverdeling  in  de  humane  huidlagen  in  volume  procenten  is  circa  2%:  4%:  94% 
(11/560:  22/560  en  527/560),  of  er  nog  metabole  activiteit  in  de  lagen  is,  is  nog  moeilijk 
vast  te  stellen.  Er  zal  in  de  komende  tijd  onderzocht  worden  wat  bekend  is  over  de 
penetratiesnelheid  van  VX  door  de  huid. 

4*2  Intraveneuze  simulaties  in  de  behaarde  cavia 

De  inputparameters  zeals  beschreven  in  voorgaande  paragrafen  worden  gebruikt  in  het 
model.  Een  listing  van  het  programma  wordt  ditmaal  niet  meegeleverd.  Het  resultaat 
van  een  van  de  eerste  simulaties  van  de  intraveneuze  toxieokinetiek  van  VX  in  de  cavia 
wordt  in  figuur  3  weergegeven. 


normal'  guinea  pig  i.v.  kinetics 


Figuur  3  Simulaties  van  de  toxieokinetiek  van  VX  na  intraveneuze  toediening  van  een  hoeveeiheid  VX 
correspondereml  mei  l,  2  en  (S  LDso (9f  IS  en  54  pg/kg)  in  de  cavia. 

Zoals  te  zien  lijkt  de  simulatie  van  1  LD50  redelijk  goed  te  Titten*  met  het  resultaat  van 
het  dierexperiment.  De  simulaties  voor  de  hogere  doses  komen  niet  meer  overeen  zoals 
deze  door  zowel  door  Boone  [7]  met  ‘gewone*  behaarde  cavia* s  werd  bepaakf 

Er  is  gekeken  naar  de  lineari teit  van  de  oppervlakte  onder  toxicokinetiek-curve  ( AUC= 
Area  Under  die  Curve)  versus  de  toegediende  dosis,  zoals  te  zien  is  in  figuur  4.  De  roze 
punten  geven  de  resultaten  weer  van  dierexperimenten  en  de  blauwe  die  van  de  simulaties. 
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Linearity 


Figmir  4  Re  present  ade  van  tie  lincariteil  van  AUC  (Area  Under  the  Curve)  versus  van  do  toegediende 

dusks.  L^e  roze  pumen  represents  re  n  d  ie  rexperi  me  at  en  en  de  blauwe  die  van  do  si  mu  I  alios. 


Het  verband  tussen  AUC  en  de  dosis  zoals  gemeteri  in.  de  dierexpenmenien  staat  haaks 
op  de  verwachting  dat  als  eenmaal  de  snelle  bindingsplaatsen  (  ACHE  en  BuChE)  verbmikt 
zijn.  de  snelheid  alleen  nog  maar  wordt  bepaald  door  de  langzame  bindingsplaatsen 
(CaE)  en  hydrolyse.  Dit  gecombineerd  zou  als  een  bi-fasisch  verloop  te  zien  moeten 
geven. 

Dat  de  werkelijkheid  toch  gecompliceerder  in  elkaar  zit  moge  duidelijk  zijn.  In  principe 
is  de  ontwikkeling  van  de  AUC  versus  de  Dosis  een  lineaire  verband  welke  het  heste 
met  een  overal  gelijk  blijvende  snelheidseonstame  beschreven  zou  moeten/kunnen  worden. 
Dit  aspect  zal  nader  onderzocht  moeten  worden. 

Een  mogelijkheid  om  hiermee  om  te  gaan  is  om  een  groter  depot  aan  a-speeifteke  bindings¬ 
plaatsen  aan  te  leggen  welke  vrijkomen  naarmate  er  meer  specifieke  bindingsplaatsen 
worden  gebruikt.  Dat  iets  dergelijks  mogelijk  zou  plaatsvinden  is  een  aantal  jaren 
bekeken  door  aan  een  weefselhomogenaat  van  de  cavia  een  hoeveelheid  radioactief 
gelabeld  organofosfaat  (soman)  toe  te  voegen  welke  meer  was  dan  de  bekende  totale 
hoeveelheid  van  AChE,  BuChE  en  CaE  te  samenj  10].  Na  incubatie  bleek  dat  dit 
homogenaat  veel  meer  organofosfaat  kon  opnemen.  De  gedachte  was  onder  andere  dat 
door  bet  homogeniseren  er  mogelijk  meer  bindingsplaatsen  werden  vrijgemaakt  Het 
lijkt  er  nu  op  dat  het  evenwicht  dat  in  compartimenten  (lees  organen)  heerst  tussen 
beschikbare  actieve  bindingsplaatsen  en  het  depot  aan  niet  direct  besehikbare  actieve 
bindingsplaatsen  wordt  beinvloed  door  de  afname  van  deze  bindingsplaatsen  door  de 
readies  met  VX.  Om  dit  te  testen  is  een  drietal  simulates  uitgevoerd  met  dezelfde 
doses  als  hiervoor,  L  2  en  6  LD50  (9t  18  en  54  jtg/kg).  waar  voor  elke  dosis  een 
toenemende  hoeveelheid  CaE  beschikbaar  is,  zie  00k  label  4. 


Tabel  4  Variaiies  in  beschikbare  eoricen trades  CaE  (nM)  in  de  compartimenten  in  relatie  met  de  dosis 
VX. 


Compartimenten 

A  fieri  eel 

Veneus 

longen 

Arm- 

Rijk- 

Lever 

Huid 

Bleed 

Bleed 

doorbloedpp 

doorbloedfp 

ILD50  (9  ng/kg) 

440 

440 

170 

5 

170 

20.000 

5 

2LDso  (18 

650 

650 

275 

5 

2.000 

60.000 

5 

6LD50  (54  pg/kg) 

650 

650 

275 

50 

2.000 

60.000 

50 
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Een  voorbeeld  van  de  invloed  van  deze  modificatie  of  de  simulatie  van  de  intraveneuze 
toxicoklineliek  van  VX  in  de  cavia  wordt  in  figuur  5  weergegeven. 


'normal'  guinea  pig  Lv,  kinetics 


Time  (min) 


Figuur  5  Si  mulattos  van  de  toxicokinetiek  van  VX  na  irnra  veneuze  toediening  van  een  hoevedheid  VX 
vergelijkbaar  mel  ],2en6  LDso  (9,  18  en  54  pg/kg)  in  de  cavia,  hierbij  nemen  de  concentralies 
CaH  toe  met  de  dosis  VX, 


Dit  concept  van  ‘automatisch'  toenemende  beschikbare  boeveelheden  CaE  voor 
ontgifting  zal  in  bet  huidige  model  nog  niet  worden  meegenomen. 

4.3  Intraveneuze  simulates  in  de  haarloze  cavia 

Door  Van  der  Sehans  et  ak  [2]  is  de  toxicokinetiek  van  VX  in  de  haarloze  cavia’ s 
onderzocht  voor  zowel  intraveneuze  toediening  alsook  de  percutane  toediening.  Een 
van  de  zaken  die  opvalt  is  het  verschil  in  de  letale  dosis  (i.v.)  welke  voor  de  haarloze 
cavia  werd  gevonden:  28  |ug/kg  versus  9  ftg/kg  voor  die  in  de  behaarde  cavia.  Dit  verschil 
kan  mogetijk  worden  veroorzaakt  door  een  hoger  aantal  bindingsplaatsen  in  de  haarloze 
cavia  ten  opzichte  van  de  behaarde  cavia.  In  hoeverre  deze  verhouding  wordt  doorgevoerd 
in  zowel  de  snelle  als  de  langzame  bindingsplaatsen  is  onbekend. 

Uitgaande  van  eenzelfde  hoevedheid  bindingsplaatsen  voor  dit  dier  als  die  in  oorspronkelijk 
in  de  behaarde  cavia  is  gevonden/gebruikt  levert  een  simulatie  op.  die  zelf  nog  slechter 
Tit*  dan  in  de  behaarde  cavia  (resultant  niet  getoond).  Indien  ook  voor  dit  dier  de 
aantal len  bindingsplaatsen  worden  opgehoogd  zoals  die  zijn  gebruikt  voor  2  LDso  in  de 
behaarde  cavia  (label  4)  dan  lijken  deze  simulates  redelijk  overeen  te  komen  met  de 
gemeten  toxicokinetiek  data,  zie  ook  figuur  6. 
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Figuur  6  Simulaties  van  de  toxicokinctiek  van  VX  na  intraveneuze  toedicning  van  een  hueveclheid  VX 
conespondereod  met  l  en  2  IDs>  (28  en  54  ng/kg)  in  do  haarloze  eavia.  hietbij  zijn  de  opgehoogde 
concent ralies  Cali  toegepast. 


Het  verschil  in  de  relatie  tussen  de  AUC  en  de  toegediende  dosis  lijkt  door  deze 
verhoging  van  CaE  concentrates  in  de  diverse  compartimenten  kleiner  te  worden  in  de 
haarloze  cavia.  tevens  ziet  emaar  uit  dat  een  redelijk  lineaire  respons  wordt  verkregen. 


I'iguur  7  Rcprcsentaiie  van  de  lincariteit  van  AUC  ( Area  Under  Lhe  Curve)  versus  van  die  de  toegediende 
dosis.  De  blauwc  punten  reprcsemeren  dierexperinienten  en  de  ro/e  punten  die  van  de  simulaties. 
hierbij  zijn  de  opgehoogde  concentraties  Call  toegepast. 


Uiteraard  is  dat  geen  'must\  het  geeft  alleen  aan  dat  het  verloop  boven  de  Letale  dosis 
(28  pg/kg)  lineair  blijkt  te  lopen.  Bij  subdetale  hoeveelheden  is  de  invloed  van  de  snelle 
bindingsplaatsen  redelijk  groot  op  het  verloop  van  de  concentratiecurve  en  daamiee  ook 
op  de  AUC.  De  lijn  loopt  niet  door  het  punt  0,  de  verwachting  is  dat  hier  eveneens  een 
bi-fasisch  verloop  te  zien  zal  zijn. 
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44  Orient  erende  percutane  simulate  in  de  haarloze  cavia 

Bij  de  percutane  experimenten  zijn  haarloze  cavia1  s  hesmet  op  de  bulk  met  een  kleine 
hoeveelheid  VX  opgelost  in  isopropanok  Onduidelijk  is  nog  of  dit  oplosmiddel  een  ml 
speelt  in  de  absorptie  van  VX  in  de  huid* 

De  opname  van  een  stof  (druppel)  door  en  in  de  huid  is  afhankelijk  van  het  expositie- 
oppervlak,  de  combi natie  van  grootte  van  de  druppel  en  de  concentrate  in  de  druppel  en 
de  penetrates  nelheid,  zie  ook  formule  5. 


KxSA 

— - x 

1000 


Celtic  ^ 


Absorbed ^ 


(5) 


SA 

Kp 

Cone  vx 
Vliq 

Absorbedvx 


-  Surface  Area  (nr) 

~  Penetratiesnelheid  (m/s) 

=  Concentratie  VX  in  de  druppel  (mol/L) 

=  druppelgrootte  (L) 

=  Qpgenomen  hoeveelheid  VX  (mol)  ->  deze  wordt  steeds  uitgerekend 

door  het  model 


Tabel  5  tnputparameters  voor  VX  absorptie  door  do  huid. 


Kp  (m/s) 

4.8*10  4 

(2,88  mm/min) 

SA(mJ) 

2*1 0  s 

(2  cm2) 

Concentratie  VX 

0.046  mol/L 

(12,3  am 

(mol/L) 

Vliq  (L) 

5*1  O'6 

(5  nl) 

Van  de  bovenstaande  parameters  zijn  or  tiric  experimented  bepaald  on  vast,  volume  van  hot  druppel tje,  do 
concent  ratio  VX  daarin  on  het  expositie  oppervlak.  De  penetratiesnelheid  is  afgeteid  nil  do  simulation 


Gebruikmakend  van  bovenstaande  mputparameters  en  die  van  de  tabellen  hiervoor 
wordt  een  simulate  verkregen  die  in  figuur  8  staat  weergegeven. 


hairless'  guinea  pig  i.v.  kinetics 


Figuur  K  Simulatie  (geirokken  lijn)  van  de  percutane  toxic okiitetiek  van  VX  bij  l  U>so  (125  fig/kg}  en 
do  meetwaarden  van  do  dierexperi menren . 
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Zo  op  hot  eersle  ge/icht  lijkt  de  simulate  reddijk  geslaugd,  Echtcr  hi j  de 
dierex  peri  men  ten  Meek  a)  dat  waarschijnlijk  (door  praktische  redencn)  le  vroeg  was 
gcsiopl  met  het  melon  van  dc  toxicokinetick  (7  hr).  Nu  is  hot  nog  onduidelijk  ofde 
curve  over  /.'n  maximum  is  en  de  concentrate  VX  in  hct  Hoed  gaat  minderen  of  dat  hot 
maximum  nog  niet  is  hereikt.  De  spreiding  in  de  met  ingen  geven  nog  gecn  eenduidig 
aniwoord. 
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5  Discussie  en  conclusie 


Ecn  behoorli jk  dee  I  van  de  benodigde  inpulparameters  is  verzameld  en  in  bet  model 
gcbruikt.  DU  PBPK  model  lijkt  na  gebruik  van  de/e  parameters  niet  direct  de  kinctick 
goed  ic  kunnen  voorspcllcn.  Na  aanpassing  van  de  lotale  hoevcelhcid  bindingsplnatsen 
oftewel  de  heschikbare  hoeveelheid  CaE  in  hepaalde  eompariimeiuen  lijkt  dil  model  de 
kinctick  van  VX  redelijk  lot  good  te  voorspclen.  De/e  heschikbare  hoevcelhcid  CaE 
lijkt  afhankelijk  tc  zijn  van  dc  aangehoden  hoevcelhcid  VX  en  aangesproken  te  worden 
wanneer  nodig.  Preliminairc  experimented  niet  tiilgevocrd  in  hel  kadcr  van  dil  project, 
li j ken  daar  ook  up  tc  wijzen. 

De  iranslormatie  van  het  PBPK  model  naareen  PBPK/PD  model  waarbij  niet  alleen  de 
toxieokinetiek  kan  worden  voorspeld  maar  tevens  het  verloop  in  de  concent  rat  ies/activite  it 
van  dc  enzymen  ACliE,  BuChE  en  CaE  kan  door  de/e  aanpak  of  gedachte  wellicht  posit  icl' 
uilvallen. 

De  voorspell  ingen  van  bet  verloop  van  de/e  enzymen  bij  suh-Ielale  doses  geven  cen  te 
snelle  en  drastische  afname  te  zien  die  mogelijk  niet  overcenkomen  met  dc  werkelijkhekk 
A  Is  vervolg  op  dil  project  in  het  kader  van  V50  2  /a]  bier  nadcr  op  in  worden  gegaan. 
Wet  lie  hi  kan  door  hel  aanbrengen  van  sub-compartimentcn  in  helangrijke  larger  organen 
uts  dc  hersenen  en  de  middenrifsptcr  een  model  worden  opgezct  dal  cen  belerc 
voorspellendc  waarde  heeft  op  het  verloop  van  AChE  en  BuChE. 

Dc  ecrsic  simulalics  van  cen  percutane  blooistclling  met  VX  bij  de  haarlo/e  eavia  waarbij 
gebruik  wordl  gemaakt  van  de  opgewaardeerde  hoevcelhcden  CaE  in  a  lie  compariimemcn 
lijkt  al  een  goede  Til'  te  geven,  Indien  dc  simutatie  lunger  wordl  gevolgd  wordt  duidelijk 
dal  het  maximum  van  VX  in  dc  hlocdconcentratie  nog  lang  niet  is  hcrcikl  (> 1 3  hr).  Dil 
onderschrijlt  het  huidige  therapie  scenario  wat  bij  cen  VX  huid  intoxicatic  gebruikelijk 
is,  namclijk  toed  ten  ing  van  injecties  met  behulp  van  auto-injectorcn  waarin  under  andere 
atropine  :/it  op  regclmaiigc  tijdsindiees  (elke  15  minuten)  na  het  optreden  van  symptomen 
lot  een  maximum  van  3.  Pas  op  de  hulppost  wordt  hel  ox  ini  rnel  lang/amc  tussenpozen 
loegediend, 

Conereel  zullen  in  V502  ecn  aantal  zaken  nader  worden  ondcr/.ochl;  de  mogeli  jkheid 
van  sub-eompartimenten  hinnen  ecn  aantal  beiangrijke  compart  imenten,  het  inbrengen 
van  cen  aantal  specifleke  toxicodynamische  compart  imenten  (hersenen  en  mogeli  jk 
midden rif)  voor  hot  uitle/cn  van  effect  en  wellicht  dc  integral ie  van  cen  therapie  model 
(alleen  voor  specifieke  compartimenten). 
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